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The near-threshold behavior of the pp invariant mass spectrum from the t/>(3686) — Y'ypp decay 
reported recently by the BESIII Collaboration is analyzed. The enhancement in the pp invariant 
mass spectrum near threshold is nicely reproduced by the pp final-state interaction based on the 
isospin averaged 1 So partial- wave amplitude as predicted by the Jiilich nucleon-antinucleon model. 
Contributions from the /2(1910) or /2(1950) mesons, as promoted in earlier works, are not needed. 



PACS numbers: 12.39.Pn; 13.25.Gv; 13.75.Cs; 25.43.+t 

Recently, the BESIII Collaboration presented data 
with improved statistics on the pp invariant mass spec- 
trum for the reaction J/-0— s>7pp, but also a first high- 
statistics measurement of the ip(3686)— >jpp decay 
The new J/ip^-^pp measurement confirmed the spec- 
tacular near-threshold enhancement in the pp invariant 
mass, found in an earlier experiment by the BES Collabo- 
ration 2] , which has been seen as evidence for a pp bound 
state or baryonium , for exotic glueball states @, [1] , 
but also simply as manifestation of the final-state inter- 
action (FSI) between the outgoing proton and antiproton 
A significant near-threshold enhancement of the 
pp invariant mass was seen also in the ^(3686) decay, 
although less pronounced than in the J/tp case. 

The BESIII Collaboration themselves interpreted their 
J/tp decay data in terms of the pp FSI proposed by us 
Q, but folded with a Breit-Wigner type resonance at 
around 1835 MeV presuming that the structure at the 
pp threshold might be related to the X(1835) resonance 
that had been observed in the reaction J/tp— >77r + 7r~r/ 
[iH [l~7l ]. see also the comments in Ref. This ob- 

ject is called X(pp) in Ref. [13] ■ For the description 
of the tp(3686) decay the same X(pp) amplitude is used 
but sizeable additional contributions from the /2(1910) 
resonance had to be invoked. Contributions of a tensor 
meson, but in this case of the /2(1950), were also advo- 
cated in the work of the CLEO collaboration 18], which 
had published data on the reaction ?/>(3686)— >jpp decay 
prior to BESIII, though with lower statistics. In both 
cases an isoscalar meson with a larger mass (Jo (2100) 
and /2(2150), respectively) has been added to explain 
the pp spectrum at higher invariant masses. 

In this report we take a closer look at those 
-0(3686)— >jpp data from the BESIII Collaboration. 
Specifically, we provide an alternative interpretation of 
the near-threshold enhancement solely in terms of the pp 
FSI, i.e. without resorting to any resonance contributions 
like the / 2 (1910) or / 2 (1950) (or the X(1835)), based on 
the very same NN interaction used by us previously in 
the explanation of the enhancement in the J/tp decay Q . 

Conservation laws for parity, charge-conjugation and 
total angular momentum severely restrict the partial 



waves in the pp system Q for such decay processes. 
Specifically, the partial-wave analysis for the J/tp de- 
cay performed in [l| suggests that the near-threshold en- 
hancement is dominantly in the J PC = h state, which 
means that the pp system should be in the partial 
wave (we use here the standard nomenclature ( 2S+1 *>Lj 
where S is the total spin and L the orbital angular mo- 
mentum). However, since the decay of the J/tp and 
■0(3686) to the jpp system involves electromagnetic pro- 
cesses, isospin is not conserved so that, in principle, any 
combination of the isospin 1 = and 1=1 components 
is allowed. Indeed, while the pp invariant mass for J/tp 
decay can be understood in terms of the FSI generated 
by the isospin 1=1 component of the NN amplitude in 
the 1 Sq state alone - at least in our work Q - the 1=1 
and 1 = channels can occur with different weights in 
case of the ^(3686) decay. 

The ^(3686)— >^pp decay rate is given by 
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where the Kallen function A is defined by A(x, y, z) = 
((x — y — z) 2 — Ayz)/Ax , M = M p p is the invariant mass 
of the pp system, fl p is the proton angle in that system, 
while fi 7 is the 7 angle in the -0(3686) rest frame. After 
averaging over the spin states and integrating over the 
angles, the differential decay rate is 
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The quantity A in Eqs. ([T]) and @ stands for the total 
■0(3686)— >7pp reaction amplitude and is dimensionless. 

We assume again the validity of the Watson-Migdal 
[T9I I20I ] approach for the treatment of the FSI effect. It 
suggests that the reaction amplitude for a production 
and/or decay reaction that is of short-ranged nature can 
be factorized in terms of an elementary (basically con- 
stant) production amplitude Aq and the pp scattering 
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FIG. 1: The pp mass spectrum from the decay 0(3686) — >"/pp. 
The circles show experimental results of the BES Collabora- 
tion [lj] . The solid line is a calculation using the pp scattering 
amplitude squared (\(Ti = i + T I=0 )/2\ 2 ) predicted by the NN 
model A(OBE) 21, 22] for the 1 So partial wave, appropriately 
normalized to the data, cf. Eq. ([3)l . The dashed (dash-dotted) 
curve are corresponding results using the I = 1 (J = 0) am- 
plitudes alone while the dotted line is the spectrum obtained 
from Eq. © by assuming a constant reaction amplitude A. 
Those curves were normalized so that they all coincide at 
M DS - 2m„ » 60 MeV. 
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FIG. 2: The pp mass spectrum from the decay 0(3686) — yypp. 
The circles show experimental results of the CLEO Collabo- 
ration [lj| . For notation of curves, see Fig. [1] 

amplitude T of the particles in the final state so that 

A(M PP ) kAq-N- T(M pP ) , (3) 

where N is an arbitrary normalization factor, see e.g. 
Ref. Q for further details. As in our investigation of 
the J/-0 decay we employ the amplitudes pred icted by 
the NN model A(OBE) published in Refs. [H [H[, and 
we assume that the FSI effects in the f/ ; (3686) decay are 
likewise dominated by the partial wave. 

Our results are presented in Fig. [T] together with the 
data from Ref. Q. As expected from the curves pre- 



sented in Fig. 3 of Ref. pj, the NN amplitude in the 
1=1 channel, which successfully describes the rather 
strong enhancement detected in the reaction J/0— >7j>p 
0, overestimates the energy dependence seen in the 
0(3686) case, cf. the dashed curve. On the other 
hand, the result based on the isospin averaged amplitude, 
(Tj=i + T I=Q )/2 = T pP , shown in Fig. Q]by the solid line, 
agrees rather nicely with the energy dependence found 
in the experiment. With an appropriately chosen nor- 
malization, cf. Eq. Q, the data are well reproduced 
from the pp threshold up to excess energies of about 150 
MeV. In particular, the \ 2 1S 8-7 for the 15 data points 
shown in Fig. Q] while it is 22.2, i.e. more than twice as 
large, for the pure 1=1 amplitude. We also include the 
result based on FSI effects due to the pp amplitude in 
the 1 = channel alone (dash-dotted curve) and we indi- 
cate the pure phase-space behaviour by the dotted curve. 
The latter is obtained by using a constant amplitude A 
in Eq. ^j. The \ 2 value for the pure 1 = amplitude is 
11.9. The one for the phase-space curve amounts to 60 
which is a clear indication that the measured invariant 
mass spectrum does not exhibit a phase-space behaviour 
near threshold. All those curves are normalized to the 
solid curve at M pp — 2m p rs 60 MeV in order to facilitate 
a comparison of the differences in the energy dependence. 

Based on those findings we do not see any need here 
to invoke further more substantial contributions coming 
from any .f 2 (1910) or / 2 (1950) mesons, say, as done in 
Refs. p], [Hj], in order to explain the data. 

In Fig. [2] our results are compared with the data 
obtained by the CLEO Collaboration p8|. Given the 
poorer statistics it is difficult to infer the actual en- 
ergy dependence of the pp invariant mass. Still, the 
solid curve corresponding to the FSI effects based on 
Tpp = (Tj = i + Tj = o)/2, again appropriately normalized, 
comes closest to the trend exhibited by the data. But 
in terms of the x 2 the differences are marginal. All re- 
sults including FSI effects yield x 2 /data « 20/30 while 
the pure phase space amounts to 30/30. 

Finally, in Fig. [3] the results for BESIII are displayed 
again, however this time in terms of the modulus squared 
of the amplitude A. Here the curves correspond directly 
to the (appropriately normalized) scattering amplitude 
squared (|T| 2 ) predicted by the NN model A(OBE) [HI, 
|22j for the 1 S'o partial wave. The symbols indicate the 
experimental values of \A\ 2 , obtained from the BESIII 
data p| via dividing the latter by the kinematical factors 
according to Eq. (0). 

In summary, we have analyzed the near-threshold 
behavior of the pp invariant mass spectrum from the 
-0(3686)— >7pp decay reported recently by the BESIII 
Collaboration within the Watson-Migdal approach. Al- 
though in this reaction there is definitely an enhancement 
in the near-threshold region as compared to the phase- 
space behavior, it is much less pronounced than what was 
found for the corresponding reaction J / ip^-^/pp. The en- 
hancement is nicely reproduced by the pp final-state in- 
teraction based on the isospin averaged 1 Sq partial- wave 
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FIG. 3: Invariant 0(3686) — yypp amplitude \A\ 2 as a function 
of the pp mass. The circles symbolize the experimental values 
of | ,4 1 2 extracted from the BES data via Eq. ©. The 
curves are the appropriately normalized scattering amplitude 
squared, |T| 2 , predicted by the NN model A(OBE) [U 
for the So partial wave. For notation of curves, see Fig. [T] 

amplitude as given by the Jiilich NN model. In par- 



ticular, any more substantial contributions from tensor 
mesons li ke / 2(1910) or /2(1950), as advocated in earlier 
works [H [13, are not required. 

Note that we have used here the same NN amplitudes 
as in our study of the J/ip decay [9(. In the Jji\> case 
the FSI provided by the 1=1 component alone led to an 
agreement with the measured near-threshold pp invariant 
mass spectrum. Clearly, the mechanisms for the decay of 
the J/ij) and "0(3686) mesons into jpp should be differ- 
ent so that different admixtures of the two isospin com- 
ponents in the final pp state have to be expected. Only 
dedicated microscopic calculations, which hopefully will 
be performed in the future, can allow to shed light on the 
details of the reaction mechanisms. 
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